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Abstract. We report 250 GHz (1.2 mm) observations of redshift > S.8 quasars from the Palomar Sky Survey 
(PSS) sample, using the Max-Planck Millimetre Bolometer (MAMBO) array at the IRAM SO-metre telescope. 
Eighteen sources were detected and upper limits were obtained for 44 with Sa flux density limits in the range 
1.5-4 mjy. Adopting typical dust temperatures of 40-50 K, we derive dust masses of a few 10^ Mq and far-infrared 
luminosities of order 10^^ Lq. We suggest that a substantial fraction of this luminosity arises from young stars, 
implying star formation rates approaching 10 Moyr"^ or more. The high millimetre detection rate supports 
current views on a connection between AGN and star forming activity, suggesting a parallel evolution of the 
central black hole and of the stellar core of a galaxy, although their growth-rate ratio seems higher than the 
mass ratio observed in nearby galaxies. The observed, exceptionally bright objects may trace the peaks of the 
primordial density field, the cores of future giant ellipticals. 

Key words. Galaxies: formation - Galaxies: starburst - Galaxies: high-redshift - Quasars: dust emission - 
Cosmology: observations - Submillimeter 



'H 1. Introduction blank field surveys at 850 ^m with SCUBA at the JCMT, 

and at 1.2 mm with MAMBO at the IRAM 30-m tele- 

^ : Understanding the history of the formation of stars and ^^^^^ j^^ve by now revealed over one hundred faint (2 to 

. - . massive black holes during the early stages of galaxy evo- ^q mJy) background sources, the majority of which re- 

lution is one of the great challenges in our quest to under- ^^^^^ invisible in deep optical and near-IR images. This 

stand the history of our Universe. Recent COBE, SCUBA, ^^^^^ conventional, spectroscopic redshift determinations 

and MAMBO observations of the far-infrared to miUime- ^f ^his population mostly impossible, 
tre extragalactic background show that most of the energy 

generated in star formation at high redshift is absorbed Global estimates, based on photometric redshifts in 
by dust, which re-emits the energy at far-infrared (far-IR) the visible/IR (e.g., Ivison et al. 2000a; Small et al. 2001) 
wavelengths, a spectral range that is red-shifted into trans- qj. (sub)niiUimetre/radio (Carilh & Yun 1999) , on the far- 
parent submiUimetre and millimetre atmospheric windows IR/submiUimetre background intensity and on theoretical 
(e.g., Gispert et al. 2000; Ivison et al. 2000a; Blain 2001; modelHng (Dole et al. 2000 and references therein), show 
Bertoldi et al. 2000a,b and Carilh et al. 2001a). that the peak of star formation in the Early Universe 
In the early stages of formation of galaxies, a large probably occured in the range z ^ 1 — 3, if we assume 
fraction of all stars apparently formed in strong bursts that much of the far-IR luminosity arises from star forma- 
within highly obscured regions. At optical and near-IR tion. It is unclear though, to what extent nucleus activ- 
wavelengths, much of the star formation activity in such ity contributes to the luminosity of high-redshift starburst 
objects remains not directly visible. Thus from the ground, galaxies, and thereby to the far-IR background. The two 
only (sub)millimetre observations provide a comprehen- brightest objects found in SCUBA and MAMBO blank 
sive measure of the energy generated in such objects. Deep field surveys (Knudsen et al. 2001; Bertoldi et al. 2000b) 
are in fact quasars, and two of the three SCUBA sources 
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signatures (Lilly et al. 1999; Smail et al. 1999; Ivison ct 
al. 2000a) 

To understand the relation between the formation of 
black holes, and the formation of galaxies and their bulges, 
it is necessary to study (aside from the (sub)millimetre 
background sources) the emission properties of distant ob- 
jects identified at X-ray, optical, near- and mid-infrared, 
and radio wavelengths. At very high redshifts (z > 4), 
where no blank field millimetre/submillimetre source has 
yet been spectroscopically identified, the relation between 
star formation and AGN is best studied through targeted 
(sub)millimetre observations of optically selected lumi- 
nous QSOs and of radio galaxies. It is still unclear to what 
extent the thermal emission of radio-quiet QSOs is pow- 
ered by starbursts, or by the black hole accretion. Just as 
for local ULIRGs and Seyfert galaxies (Genzel et al. 1998), 
there is increasing evidence that in high-redshift AGN a 
substantial part of the thermal emission comes from star- 
bursts. 

Well before any SCUBA or MAMBO deep imaging 
surveys, we were able to detect seven powerful high- 
redshift millimetre sources through targeted observations 
of ultra-luminous z > 4 quasars using the 7- and 19- 
element bolometer arrays at the IRAM 30-metre telescope 
(McMahon et al. 1994; Omont et al. 1996a) and the JCMT 
(Isaak et al. 1994). Evidence for the starbust origin of their 
millimetre emission arose from the determination of their 
spectral energy distribution, including 850-450 /im obser- 
vations at the JCMT (Isaak et al. 1994, McMahon et al. 
1999) and 350 /im measurements at the CSO (Benford et 
al. 1999), and from the detection of CO emission in three 
objects (Ohta at al. 1996; Omont et al. 1996b; Guilloteau 
et al. 1997, 1999). 

We here report on a substantial extension of the sur- 
veys for 1.2 mm continuum emission from z > 4 ultra- 
luminous quasars with 18 detections, which multiplies by 
three the number of detections. Results of 850 /im obser- 
vations of a part of the present sample are presented in 
Isaak et al. (2001). 

2. Observations and Results 

The observations were made in February and March 2000 
with the 30- metre IRAM telescope at Pico Veleta (Spain), 
using the Max-Planck Millimetre Bolometer (MAMBO; 
Kreysa et al. 1999). MAMBO is a 37-element bolome- 
ter array operating at 300 mK. It is sensitive between 190 
and 315 GHz, with half-power sensitivity limits at 210 and 
290 GHz, and an effective bandwidth center for steep spec- 
tra at 250 GHz. The bolometer feed horns are matched to 
the telescope FHWM beam at 1.2 mm of 10'.'6, and they 
are arranged in an hexagonal pattern with a beam sepa- 
ration of 22". The sources were observed with the array's 
central channel, using the standard on-off observing mode 
with the telescope secondary chopping in azimuth by 50" 
at a rate of 2 Hz. The pointing was monitored regularly on 
nearby sources and was found to be stable typically within 
~ 2". Gain calibration was performed using observations 



of the planets Mars and Uranus as well as the asteroid 
Ceres. A calibration factor of 12500 counts per Jansky is 
adopted, which we estimate to be accurate to within 20%. 
The data were reduced using the MOPSI software pack- 
age (Zylka 1998). The point-source sensitivity of MAMBO 
at 1.2 mm is 20 to 40 mJy/\/Hz, depending on weather 
conditions and the effectiveness of sky-noise subtraction. 
In order to remove the strong sky signal from the un- 
derlying source signal, we subtract from each channel the 
weighted, average correlated signal from its surrounding 
channels. Subtracting correlated noise has shown to effec- 
tively reduce the sky-noise by typical factors of 2 to 3. 
Because the channels are separated on the array by two 
beams, point sources do not produce correlated signals be- 
tween neighboring channels. Since the source signal is not 
affected by the sky-noise subtraction (which is confirmed 
by extensive tests), the error on the final source signal is 
given by the remaining noise level of the source channel. 
Systematic errors however are introduced by bad point- 
ing, focus, and gain variations that lead to uncertainties 
in the flux calibration. We have eliminated observations 
where we judge such errors to be significantly larger than 
the errors due to the normal system- and sky-noise. The 
remaining flux calibration uncertainty are of order 10%. 

Most of our sources were observed several times on dif- 
ferent dates, and we were able to check for consistency in 
the resulting fluxes. We noticed that several observing pe- 
riods yielded source fluxes which were altogether inconsis- 
tent with observations of the same objects on other dates. 
We suspect that during such periods, either the telescope 
pointing, optics, or an unstable atmosphere lead to a sig- 
niflcant sensitivity loss, and we therefore ignored observa- 
tions during such suspect periods. For a reliable flux de- 
termination, a source should be observed at least on two 
different dates and yield consistent results, or a reliable, 
nearby flux calibrator should have been observed along 
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Fig. 1. Redshift distribution for the PSS quasars which 
were detected at 1.2 mm (filled histogram) and for those 
not detected (dotted line). 
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Table 1. Quasars detected at 1.2 mm. 



Source 


2 


Mb^ 


R.A. 


Dec. 


Flux Density 


Quality* 


Time 




Comments" 








(J2000.0) 


(mjy, ±la) 




(sec) 






PSS J0209+0517 


4.18 


-28.1 


02 09 44.7 


05 17 13.3 


3.3±0.6 


A 


2240 






PSS J0439-0207 


4.40 


-27.5 


04 39 23.1 


-02 07 01.7 


2.3±0.7 


A 


1940 


RL^ 




PSS J0808+5215 


4.44 


-28.7 


08 08 49.4 


52 15 15.2 


6.6±0.6 


A 


1420 


RQ' 


smm 


PSS J1048+4407 


4.40 


-27.4 


10 48 46.6 


44 07 12.7 


4.6±0.4 


A 


3100 


(b) (c) W As Comp smm 


PSS J1057+4555 


4.12 


-28.8 


10 57 56.3 


45 55 53.3 


4.9±0.7 


A 


1930 


(c) W As Comp smm 


BRB1117-1329 


3.96 


-28.2 


11 20 10.2 


-13 46 26.2 


4.1±0.7 


A 


1420 


(a) (b) BAL smml 


BR B1144-0723 


4.15 


-27.7 


11 46 35.6 


-07 40 05.5 


6.0±0.7 


A 


2000 


(a) (b) BAL smml 


PSS J1226+0950 


4.34 


-27.4 


12 26 23.8 


09 50 04.7 


2.8±0.7 


A 


1410 






PSS J1248+3110 


4.35 


-27.6 


12 48 20.2 


31 10 44.3 


6.3±0.8 


A 


2050 


RQ' 


smm 


PSS J1253-0228 


4.00 


-27.2 


12 53 36.4 


-02 28 08.3 


5.5±0.8 


A 


1290 


RQ' 


Comp 


PSS J1317+3531 


4.36 


-27.5 


13 17 43.1 


35 31 32.3 


3.7±1.1 


B 


1160 


RQ' 


RQ^ W As 


PSS J1347+4956 


4.56 


-28.3 


13 47 43.3 


49 56 21.3 


5.7±0.7 


A 


1650 






PSS J1403+4126 


3.85 


-26.8 


14 03 55.7 


41 26 16.2 


1.5±0.5 


A 


2300 






PSS J1418+4449 


4.32 


-28.6 


14 18 31.7 


44 49 37.6 


6.3±0.7 


A 


1290 


RQ' 


smm 


PSS J1535+2943 


3.99 


-27.1 


15 35 53.9 


29 43 13.0 


1.9±0.6 


A 


1410 






PSS J1554+1835 


3.99 


-26.6 


15 54 09.9 


18 35 51.0 


6.7±1.1 


A 


1170 


RQ' 




PSS J1555+2003 


4.22 


-27.4 


15 55 02.6 


20 03 24.4 


3.1±0.6 


A 


1410 






PSS J1646+5514 


4.10 


-28.7 


16 46 56.5 


55 14 46.0 


4.6±1.5 


B 


380 


RQ' 


smm 


PSS J1745+6846 


4.13 


-27.0 


17 45 50.4 


68 46 20.5 


2.5±0.7 


B 


1530 






PSS J1802+5616 


4.18 


-26.9 


18 02 48.8 


56 16 49.6 


2.8±0.9 


A 


1540 






PSS J2322+1944 


4.11 


-28.1 


23 22 07.2 


19 44 23.0 


9.6±0.5 


A 


1930 


RQ' 


smm 



Notes: Redshifts and coordinates of observed positions for the PSS quasars are from G. Djorgovski's web page (see text). 

Redshifts and coordinates for the APM BRI quasars are from Storrie-Lombardi et al. (1996). Note that new astrometry has led 

to a redetermination of the position of PSS J1057+4555, which is now listed as 10 57 56.39, +45 55 51.97. For PSS J1048+4407, 

Isaak et al. (2001) point out that the correct position is 10 48 46.64 +44 07 10.9. 

' Absolute B-band magnitudes have been calculated as described in Isaak et al. (2001), for Hq = 50 kms^^ Mpc^^ and go=0.5. 

* Quality of the observations: A is good and B is poor. 

' (a) Omont et al. (1996a) (b) Guilloteau et al. (1999) (c) see text for radio properties. RQ=radio-quiet, RL=radio-loud: the 

subscripts 1, 2 and 3 refer to Carilli et al. (2001b), Stern et al. (2000) and Isaak et al. (2001), respectively; BAL=Broad 

Absorption Lines; W=Weak (and broad) emission lines; As= Associated absorption system; Comp=companion(s); smm and 

smml refer to 850 /im detections reported by Isaak et al. (2001) and McMahon et al. (1999), respectively. 



with the source. These considerations were used to define 
a quality factor which we list with the fluxes in Table 1. 

The high-redshift quasars were selected from the mul- 
ticolor Palomar Digital Sky Survey (DPOSS: Kcnnefick et 
al. 1995a and b; Djorgovski et al. 1999), that were avail- 
able at the time from G. Djorgovski's web page: 



http: //astro. caltech.edu^george/z4.qsos 



In order to study the most luminous objects, we arbi- 
trarily limited the observations to sources brighter than 
Mb — —27.0 with a few exceptions. Altogether, 62 PSS 
sources were observed out of the 90 PSS quasars listed by 
G. Djorgovski at the time. 

Eighteen new sources were detected at levels > 3 cr, of 
which 12 sources were detected at levels <; 4cr, as out- 
lined in Table 1. Figure |l| shows that their redshifts are 
miiformly distributed between z — 3.9 and 4.5. Figure ^ 



displays the observed 1.2 mm flux density 6*250 versus 
Mb, the optical absolute magnitude in the B rest- frame 
band, which represents well the bolometric magnitude, 
Lboi (Sec. 3.6). It shows that there is no strong correlation 
between Mb and 5250 ■ 

Nine of the quasars are detected with a 1.2 mm flux 
density ^ 5mJy. PSS J2322+1944, with a flux density of 
9.6±0.6 mJy, is the strongest millimetre source at z > 4 
just after BR1202— 0725, which has a 1.2 mm flux den- 
sity of 18.5±2 mJy (derived from Guilloteau et al. 1999). 
Both sources are exceptional and are indeed brighter 
than any other high-redshift source known, except the 
most strongly lensed sources, H1413+117 (Cloverleaf), 
FIRAS10214+4724 and APM08279+5255 (however, see 
also the case of BRI B1335-0417, Guilloteau et al. 1999). 
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Table 2. Quasars with upper limits at 1.2 mm. 



Source 


2 


Mat 


R.A. 


Dec. 


Flux Density 


Quality * 


Time 


Comments 








(J2000.0) 


(mjy, ±lcr) 




(sec) 




PSS J0014+3032 


4.40 


-28.1 


00 14 43.0 


30 32 03.3 


-0.5 ±0.9 


B 


780 




PSS J0121+0347 


4.13 


-27.6 


01 21 26.1 


03 47 07.2 


1.7±1.2 


B 


780 


RL^, RL^ 


PSS J0121+3453 


4.22 


-27.6 


01 21 11.5 


34 53 00.0 


0.7 ±1.0 


B 


780 




PSS J0131+0633 


4.43 


-27.8 


01 31 12.2 


06 33 40.0 


1.3 ± 1.2 


B 


780 




PSS J0134+3307 


4.52 


-27.4 


01 34 21.5 


33 07 55.4 


3.5 ±1.5 


A 


780 


RQ' 


PSS J0211+1107 


3.99 


-27.1 


02 11 20.1 


11 07 16.0 


-1.0 ±0.9 


A 


1420 


RL3 


PSS J0452+0355 


4.38 


-27.6 


04 52 51.5 


03 55 57.0 


0.2 ±0.8 


A 


1300 


smm 


PSS J0747+4434 


4.42 


-28.5 


07 47 49.7 


44 34 16.0 


0.1 ± 1.2 


A 


780 


RQ^ 


PSS J0852+5045 


4.20 


-28.2 


08 52 27.3 


50 45 10.9 


0.6 ±0.7 


A 


1300 




PSS J0926+3055 


4.19 


-29.5 


09 26 36.3 


30 55 04.9 


1.0 ±0.8 


A 


2330 




PSS J0950+5801 


3.97 


-28.3 


09 50 13.9 


58 01 38.0 


0.4 ±0.9 


A 


780 




PSS J0955+5940 


4.34 


-28.9 


09 55 11.3 


59 40 30.5 


-1.3 ± 1.1 


A 


900 




PSS J0957+3308 


4.25 


-28.3 


09 57 44.5 


33 08 20.4 


0.1 ± 1.2 


B 


640 




PSS J1026+3828 


4.18 


-27.8 


10 26 56.7 


38 28 44.9 


0.4 ±0.7 


B 


900 




PSS J1058+1245 


4.33 


-28.8 


10 58 58.4 


12 45 54.8 


0.1 ±0.9 


B 


1290 




PSS J1159+1337 


4.08 


-28.8 


11 59 06.6 


13 37 37.6 


0.9 ±0.8 


B 


1380 


RQ' 


PSS J1315+2924 


4.18 


-27.4 


13 15 39.8 


29 24 39.1 


0.7 ±0.7 


A 


900 




PSS J1326+0743 


4.17 


-28.4 


13 26 11.8 


07 43 57.5 


-0.3 ±0.6 


A 


1800 




PSS J1339+5154 


4.08 


-27.2 


13 39 13.0 


51 54 03.8 


1.4 ±0.6 


B 


1540 




PSS J1401+4111 


4.01 


-27.9 


14 01 32.8 


41 11 50.9 


-0.4 ± 1.2 


B 


520 




PSS J1430+2828 


4.30 


-27.1 


14 30 31.6 


28 28 34.3 


3.2 ±1.2 


B 


650 




PSS J1432+3940 


4.28 


-28.2 


14 32 24.8 


39 40 24.4 


-0.1 ±0.7 


A 


1420 




PSS J1435+3057 


4.35 


-28.5 


14 35 23.5 


30 57 22.1 


2.4 ± 1.6 


A 


390 


RQ2 W 


PSS J1443+2724 


4.42 


-27.3 


14 43 31.2 


27 24 36.9 


-0.1 ±0.5 


A 


2190 


RQ^ 


PSS J1443+5856 


4.27 


-28.5 


14 43 40.7 


58 56 53.3 


0.7 ±0.6 


B 


2320 




PSS J1500+5829 


4.22 


-27.8 


15 00 07.7 


58 29 37.7 


0.4 ±0.9 


A 


640 




PSS J1506+5220 


4.18 


-27.7 


15 06 54.6 


52 20 04.6 


3.2 ±1.2 


A 


780 




PSS J1531+4517 


4.20 


-27.4 


15 31 29.9 


45 17 07.9 


0.9 ±0.7 


B 


1940 




PSS J1543+3417 


4.39 


-29.4 


15 43 40.4 


34 17 44.5 


0.3 ±0.6 


A 


1290 




PSS J1615+1803 


4.01 


-27.6 


16 15 22.9 


18 03 56.0 


-0.1 ±0.7 


A 


1420 




PSS J1618+4125 


4.21 


-27.4 


16 18 22.7 


41 25 59.7 


1.2 ±1.1 


B 


500 


RQ' 


PSS J1633+1411 


4.35 


-27.8 


16 33 19.6 


14 11 42.3 


1.1 ±1.2 


B 


640 




PSS J1721+3256 


4.03 


-27.1 


17 21 06.7 


32 56 35.9 


-0.8 ±0.8 


A 


1160 


RQ2 Comp 


PSS J2134+0817 


4.** 


-27.1 


21 34 43.2 


08 17 28.6 


0.3 ±0.9 


A 


780 




PSS J2154+0335 


4.46 


-27.5 


21 54 06.7 


03 35 39.2 


0.4 ±0.7 


A 


1530 




PSS J2155+1358 


4.26 


-28.1 


21 55 02.2 


13 58 25.5 


0.6 ±0.8 


B 


910 




PSS J2203+1824 


4.38 


-28.7 


22 03 43.4 


18 28 13.1 


0.4 ±0.9 


B 


780 




PSS J2238+2603 


4.03 


-28.9 


22 38 41.6 


26 03 45.3 


2.5 ±1.3 


B 


750 




PSS J2241+1352 


4.46 


-28.0 


22 41 47.8 


13 52 02.0 


-0.0 ±1.3 


A 


510 




PSS J2244+1005 


4.04 


-26.9 


22 44 05.4 


10 47 38.6 


-1.5 ±1.3 


A 


520 




PSS J2256+3230 


4.04 


-26.0 


22 56 10.4 


32 30 18.5 


-1.7 ±0.8 


B 


780 


RL^ 


PSS J2315+0921 


4.52 


-28.1 


23 15 59.1 


09 21 43.6 


-0.4 ± 1.1 


B 


910 




PSS J2323+2758 


4.18 


-27.3 


23 23 40.9 


27 58 00.4 


1.7±1.1 


A 


780 




PSS J2344-f 0342 


4.30 


-28.1 


23 44 03.2 


03 42 26.4 


0.1 ±0.8 


A 


1540 





Note: Same notations as in Table 1. 
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Fig. 2. Observed 1.2 mm flux density versus Mb, the op- 
tical absolute magnitude in the B rest-frame band. QSOs 
detected at 1.2 mm (Table 1) are shown as filled squares. 
The triangles represent 2 a upper limits for sources not 
detected at 1.2 mm (Table 2) - plotted are the 2 a values 
or the flux density plus a if ^250 > 2 a. 



Interestingly, all these single-dish measurements (ex- 
cept at 850 /im for BR B1144— 072) are incompatible, for 
point sources, with the results of the IRAM interferome- 
ter measurements at 1.35 mm reported by Guilloteau et 
aL (1999). If PSS J1048+4407, BR B1117-1329 and BR 
B1144— 0723 are unresolved by the interferometer beam 
(« 3" X 2"), their 1.35 mm flux densities should be limited 
to 0.25 ± 0.68 mJy, 0.64 ± 1.47 mJy and 0.54 ± 0.78 mJy, 
respectively. The most natural explanation of the poor 
agreement between the single dish and the interferome- 
ter would be that the dust millimetre emission of these 
three sources is extended with respect to the interferom- 
eter beam. Such extensions have been reported for other 
high redshift sources (Omont et al. 1996b; Papadopoulos 
et al. 2000, 2001; Richards 2000). However, this needs con- 
firmation for PSS J1048+4407, BR B1117-1329 and BR 
B1144-0723. 

Forty-four PSS quasars were not 3 a detected at 
1.2 mm, with typical 3 a fiux density upper limit of 1.5- 
4 mJy. Table 2 lists their names, redshifts, and the 1.2 mm 
fiux densities with ±1 a errors. 

We have detected 30% of the 62 new sources observed. 
Such a detection rate is consistent with previous studies, 
in particular with that of Omont et al. (1996a) on a three 
times smaller sample of similar brightness APM BRI sur- 
vey quasars. It is also consistent with the parallel, deeper 
MAMBO study of a sample of 41 SDSS quasars (Carilh 
et al. 2000, 2001a). It should be stressed that such bright 
z > 4 quasars with M^ < —27 are rare, with an average 
surface density of one per 50-100 deg^ (Storrie-Lombardi 
et al. 2001). The number of z > 4 quasars detectable in 
the (sub) millimetre with fiux densities ^ 5mJy is thus 
less than one per 100 deg^. 



In addition, three quasars (PSS J1048-F4407, BR 
B1117-1329 and BR B1144-0723) which were previously 
detected at the 30-metrc telescope, but were not detected 
at similar levels in observations with the Plateau de Bure 
interferometer (Omont et al. 1996a, Guilloteau et al. 1999 
Table 2), were remeasured. Their new 1.2 mm flux densi- 
ties are between 4 and 5 mJy (Table 1), confirming the pre- 
vious 30-metre data (3±1, 4.09±0.81, and 5.85±1.03mJy, 
respectively, at 1.25 mm). These values are also compati- 
ble with the 850 /im fiux densities measured with SCUBA 
at the JCMT: 12 ± 3mJy for PSS J1048+4407 (Isaak et 
al. 2001), 13 ± 1 mJy and 7 ± 2 mJy for BR B1117-1329 
and BR Bl 144-0723, respectively (McMahon et al. 1999). 
This consistent set of data makes the fiux densities mea- 
sured at 1.2 mm and 850 fim with single dishes very ro- 
bust for the three sources. For BR B1144— 0723 a new 
measurement at 850 /im with a better sensitivity would 
be desirable. As a consistency check, this source was re- 
measured in February 2001 with the 30- metre telescope 
at 1.2 mm yielding a result which is perfectly consistent 
with the measurements done in 2000. The value of the flux 
density listed in Table 1 includes this recent measurement. 



3. Discussion 

3.1. The case for thermal dust emission 

There is little doubt that, for nearly all the high-redshift 
quasars detected in the millimetre, this emission is due 
to dust, and not to synchrotron radiation. This is shown 
either by their large and positive submillimetre/millimetre 
spectral index, which is characteristic of dust emission, 
or by the fact that they are not radio-loud at 1.4 GHz 
or 5 GHz. Both indicators have well established the dust 
origin of the millimetre emission of the z > 4 quasars 
detected until now at 1.2 mm (Isaak et al. 1994; Benford et 
al. 1999; McMahon et al. 1999; Carilh et al. 1999; Kawabe 
et al. 1999; Priddey and McMahon 2001). 

The characteristics of the dust emission of previously 
detected sources have been derived by combining obser- 
vations made at various (sub)millimetre wavelengths. For 
instance, by comparing 350 /im measurements of six high- 
z quasars with 1.3 mm - 450 /im fiux densities, Benford 
et al. (1999) showed that most of the dust contributing 
to the millimetre/submillimetre emission is at a temper- 
ature close to 50 K, with a power-law emissivity index 
/3==1.5±0.2. Recently, Priddey & McMahon (2001) com- 
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piled an average spectral energy distribution (SED) for 
the z > 4 quasars detected by Omont et al. (1996a) 
and derived T = 41 ± 5K and /3=1.95±0.3. Note that 
these data are not exclusive of dust at higher tempera- 
ture, as seen in H1413+117 (Cloverleaf ) , F10214+4724 
or APM08279+5255, and expected from models (e.g., 
Andreani et al. 1999). The luminosities which have been 
derived in previous studies from the SEDs of z > 4 
quasars detected in the millimetre are > 10^'^ L0 (see 
also Eq. 3) , slightly greater than the typical luminosity of 
ultra-luminous infrared galaxies such as Arp 220. 

Similarly, radio and submillimetre data begin to be 
available for a part of the broader PSS sample, especially 
for the sources that we have detected at 1.2 mm. All results 
point to dust emission. A very deep radio VLA study of 
nine PSS QSOs, among the strongest 1.2 mm sources of 
Table 1 , shows that their millimetre emission is most likely 
thermal dust emission (Carilli et al. 2001b). Six of them 
are radio-quiet (Table 1). PSS J1057+4555 is shown as a 
weak radio-loud quasar and PSS J1048-I-4407 as a weaker 
one. However, the steep radio spectrum of both sources 
makes it unlikely that their large 1.2 mm flux densities 
result from synchrotron emission. 

Stern et al. (2000) and Isaak et al. (2001) have also 
addressed the radio properties of optically selected z > 4 
quasars. From a deep 5 GHz imaging survey. Stern et al. 
have identified as radio-quiet two of the sources of Table 
1 (detected at 1.2 mm) and seven of Table 2 (undetected 
at 1.2 mm). From a comparison with the shallow NVSS 
1.4 GHz survey, Isaak et al. (2001) have identified as radio- 
loud one source detected at 1.2 mm, PSS J0439— 0207, and 
three undetected ones (Table 2). 

Seven of the 1.2 mm MAMBO PSS detections (marked 
by 'smm' in Table 1) were also detected at 850 /im by 
Isaak et al. (2001), indicating 850 /xm/1.2 mm flux ratios 
which are characteristic of dust emission, in agreement 
with their radio properties. A detailed discussion of the 
submillimetre properties of these quasars and parameters 
inferred from them is left to a later publication (Priddey 
et al., in preparation). 

3.2. Dust mass 

In the likely case that the millimetre emission is optically 
thin, the total mass of dust is directly proportional to the 
millimetre flux 



Mrf = 



5*250 D, 



(1 + z)Kd{Vr)B{Vr,Td) 



(1) 



where S'250 is the flux density at 1.2 mm (i.e., 250 GHz), Vr 
is the rest- frame frequency, Kd{i'r) is the dust absorption 
coefficient, B(i^r, Tj) is the Planck blackbody function and 
Di^ is the luminosity distance which depends on the cos- 
mological parameters (see, e.g., Eq. (1) in McMahon et al. 
1994). Hereafter, we adopt a standard Einstein-de Sitter 
cosmology with Hq — 50 kms""'^ Mpc~^ and (70=0. 5. To 
convert to a Lambda cosmology with f^M = 0.3, f^A — 0.7, 
and Hq~65 kms~^ Mpc~^, the masses and luminosities 



given here should be multiplied by « 1.4. For the median 
value of the redshift distribution of Fig. 1, z « 4.2, the 
rest-frame wavelength is R:!230^m. Equation (1) can then 
be rewritten as 

where F{Td) = {exp[hi^ / kTd] - l)/(exp[/i://fc • 50K] - 1) 
accounts for the temperature dependence. The dust absop- 
tion coefficient at 230 /im, Kd{230 fim) = 7.5cm^g~^, is the 
value extrapolated from the value at 125 /im of Hildebrand 
(1983), Kd{125 fim) ~ 18.75 cm^g^^ and a frequency de- 
pendence of /3— 1.5. Using Eq. (2) and Td ^ 50 K, the 
mass of dust of the sources detected at 1.2 mm (Table 1) 
is in the range 2 x 10^ to 10^ (7.5/Kd23o) Mq with a median 
value of - 5 X 10^ (7.5/Kd23o) Mq. 

Equation (2) is consistent with previous formulae (e.g., 
Omont et al. 1996a; Rowan- Robinson 2000). However, 
note that for Td ^ 40 — 50 K the above relation is slightly 
different from Eq. (1) in Omont et al. (1996a) where a 
dust temperature of 80 K was assumed. For T^ ^ 50 K, 
the temperature dependence of Md is roughly T^^'^- 

Note that other temperature dust components may be 
present. As quoted in Sec. 3.1, in the few quasars for which 
higher frequency measurements are available, a warm dust 
component is evident. A significant cold dust component 
could also be hidden below the spectral energy distribu- 
tion of the warm dust. In the most extreme case of dust 
at the CMB temperature at z = 4, T^ w 13 K, a given 
fiux density at 1.2 mm could arise from five times the 
amount of dust we would estimate assuming Td — 50 K. 
More realistically, we estimate that within the flux density 
uncertainties a hidden cold dust component could still be 
present with perhaps about twice the mass we estimate 
ignoring such a cold component. 

We should emphasize that the estimates of M^ depend 
on the assumed values of the dust temperature and emis- 
sivity. The latter remains highly uncertain, particularly 
for the unknown properties of dust in luminous z > 4 
quasars. Dust masses Md ~ 1 — 4 x 10^ Mq are de- 
rived by Andreani et al. (1999) in their Table 3 for the 
sources whose 1.25 mm flux values are reported by Omont 
et al. (1996a). The origin of such low values of Md is 
probably due to the high dust temperature resulting from 
AGN heating and to the adopted grain model (Granato 
& Danese 1994; Rowan-Robinson 1986) which has a large 
emissivity at millimetre wavelengths. 

3.3. Dust heating: Starburst or AGN? 

As discussed by various authors (Sanders & Mirabel 1996; 
Andreani et al. 1999; Rowan-Robinson 2000; Yun et al. 
2000; Carilli et al. 2000), for ultra luminous infrared galax- 
ies (ULIRGs) it is difficult to disentangle the dust heating 
contributions from a starburst and from the central AGN, 
even in the local universe (Genzel et al. 1998; Downes & 
Solomon 1998; Veilleux et al. 1999), where such objects are 
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bright but rare. Due to the lack of detailed spectral, geo- 
metric, and kinematic information, for sources at very high 
redshifts only indirect arguments can be made in favour 
of one or the other heating mechanism. 

An argument that favours AGN as the dominant heat- 
ing source is the large power available that could heat the 
dust even at kpc distances. Simple torus models can ac- 
count for the observed SED of AGN at small and moder- 
ate redshift (Sanders et al. 1989; Granato & Danese 1994; 
Granato et al. 1996; Andreani et al. 1999). In particu- 
lar, AGN heating seems required to explain the dust de- 
tected at higher temperature in lensed sources such as 
H1413+117, FIRAS10214+4724 and APM08279+5255, if 
not due to differential lensing (e.g., Rowan-Robinson 2000 
and references therein). 

Starbursts as a source of the far-IR luminosity are 
favoured by the large observed dust masses, ^ 10^ M©, 
which imply the presence of > 10^" Mq of gas, much of 
which may be molecular as indicated by the detection of 
CO in three z > 4 quasars (Ohta et al. 1996; Omont et al. 
1996b; Guilloteau et al. 1997, 1999). Such gas masses are 
typical for local ULIRGs (Solomon et al. 1997; Downes 
& Solomon 1998) and hyper- luminous infrared galaxies 
(Rowan- Robinson 2000). The high gas masses, the dust 
temperatures (Benford et al. 1999), and the faint radio 
emission (Yun et al. 1999; Kawabe et al. 1999; Carilli et 
al. 2000, 2001b) are consistent with starbursts with star 
formation rates, AIsf, of order 1000 MQyr~^. The total 
far-IR luminosity, LpiR, which we take to be the integrated 
luminosity at rest wavelengths > 50 /um, has been well de- 
termined from millimetre to 350 /im detections (Benford 
et al. 1999) for several objects at z > 4. 

For a modified black-body with Td — 50 K and 
P = 1.5, LpiR is related to Md through LpiR = 3.3 x 
10^2 (Md/lO^Me) (Kd23o/7.5cm2g-i) Lq - see Eq. (5) in 
McMahon et al. (1999). Therefore, from Eq. (2) 



Lfir ~ 3.6 X 10^2 (525o/mJy) Lg 



3.4. Starbursts: Implied star formation rates 



(3) 



Let us consider the case that the far-IR luminosity arises 
exclusively from starbursts, and that the stellar radiation 
is completely absorbed and reemitted by warm dust. Then 
the far-IR luminosity equals the starburst's total lumi- 
nosity. In a continuous starburst the star formation rate 
relates to the far-IR luminosity 



A^SF = Smf (Lfir/101*'Lq) Mgyr- 



(4) 



where Smf is a function of the mass function (MF), i.e., 
of the present mass composition of the stellar popula- 
tion. To determine Smf, as an example, we first consider 
stars forming with a Salpeter initial mass function (IMF), 
d log iV/d log TO = —2.35, with a low mass cutoff, to/. We 
can adopt the approximate expressions for the stellar lu- 
minosities and lifetimes given by Telesco & Gatley (1984). 
Then for a 100 Myr old continuous starburst, Smf = 1.2 




1000 



starburst age [Myr] 



Fig. 3. Star formation rate per total luminosity, Smf, 
for continuous starbursts as a function of starburst age, 
adopting different IMFs. Solid lines: Salpeter IMF with 
TO,; = O.IMq (upper solid line) and mi = 1.6Mq (lower 
solid line), with stellar lifetimes and luminosities taken 
from Schaerer et al. (1993) evolution tracks. Dashed lines: 
As for solid lines, but adopting approximate lifetimes and 
luminosities from Telesco & Gatley (1984). Dash-dotted 
line: Like upper solid line but assuming a flat IMF be- 
low 1M0. Dotted line: As upper solid line, but assuming 
a log-normal IMF: dlogiV/dlogTO = — log to. 



and 3.8, for a low mass cutoff to,; — 1.6 M© and 0.1 M©, re- 
spectively, (Fig. H). These values of Smf are lower by a fac- 
tor of two from those which are commonly adopted from 
Thronson & Telesco (1986), which were computed for a 
continuous starburst of age 2 Myr. Such a short starburst 
age is certainly inappropriate for the massive starbursts 
we observe in high-rcdshift galaxies, which more likely oc- 
cur on timescales comparable to the dynamical timescale 
of the starburst regions, lO"^ — 10* yr. 

We have alternatively computed Smf using recent stel- 
lar evolution models (Schaerer et al. 1993), considering 
a range of metallicities and IMFs. For starburst ages of 
10-100 Myr, and a more realistic flat IMF at low masses 
{d\ogN/d\ogm = -1 at to < IM©, mi = 0.1 M©), ^mf 
ranges between 0.8 and 2 (Fig. g), well in agreement with 
previous studies (see, e.g., Kennicutt 1998; Barger et al. 
2000; Rowan- Robinson 2000). Whilst the exact value of 
Smf for the z > 4 sample discussed in this paper is un- 
known, it is unlikely that it will deviate much from this 
range. 

In reality, the far-IR luminosity may differ from the 
starburst's total luminosity for two main reasons: either 
part of the stellar radiation is not absorbed by dust and 
escapes, or part of the dust heating is achieved by the 
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quasar radiation. Let Ssb be the ratio between the star- 
burst's total luminosity and Lfir- Note that Ssb is just an 
empirical definition which can depend on several parame- 
ters. It follows from Eqs. (3) and (4) that the star forma- 
tion rate of a continuous starburst relates to the 1.2 mm 
flux density, 6*250, 



AfsF ~ 360 5mf Ssb (6250 /mJy) Moyr" 



3.5. Heavy element production and FIR emission 



(5) 



The dust seen through its far-IR emission was generated 
by stars that formed during current or past star formation 
events. The mass of dust may provide an estimate of the 
minimum amount of heavy elements (i.e., > 12 amu) al- 
ready synthesized, and hence of the energy radiated from 
nucleosynthesis and of the total star formation. The en- 
ergy directly generated in the nucleosynthesis of a mass 
Mm of heavy elements is ~0.008 Mm c^. Accounting for 
the energy released in the formation of helium which is si- 
multaneously synthesized in massive stars, but not trans- 
formed in heavy elements, one may estimate that the ra- 
diative energy, Em, emitted by the stars which release a 
mass. Mm, of heavy elements into the interstellar medium, 
is of the order of 



Er, 



0.02 Mm c'' 



(6) 



(see Harwit 1999). From Eqs. (4) and (6), the factor of 
proportionality betwen the total mass of formed stars and 
of released heavy elements is 



Msf/M„ 



30 6 



MF 



(7) 



The mass of heavy elements ever released in a galaxy 
must be larger than the mass of dust we can observe 
through its FIR emission. We introduce 



Smd = Mm/Md > 1 



(8) 



as the ratio between the total mass of released heavy ele- 
ments and the visible dust mass. This ratio accounts for 
heavy elements in the gas phase and in the intergalactic 
medium, or reprocessed in new stars, and for cold dust. In 
the local interstellar medium, the mass fraction of heavy 
elements locked in dust grains is ~40%, when adopting 
a total abundance of heavy elements ^2/3 solar (Meyer 
et al. 1998) and a gas-to-dust ratio --190 (Allen 1999, 
Table 7.12). This fraction is unknown in local starburst 
galaxies, but there are no indications that it is very differ- 
ent, and the presence of a consistent amount of gas phase 
oxygen and carbon at least is attested through the detec- 
tion of CO. Adopting the Galactic dust fraction of heavy 
elements, we speculate that in high-redshift starbursts, 
Smd ~ 2.5, but that this value is uncertain by a factor of 
two in either direction. 

Some of the heavy elements produced in a starburst 
may be blown into the intergalactic medium or scattered 
throughout the galaxy to regions far from radiative energy 
sources. The dust in this component may not be heated 



sufficiently to contribute noticably to the integrated far-IR 
emission. The intergalactic ejection of enriched gas is ex- 
pected to be particularly significant in extreme starbursts 
such as those we might be observing in the high-redshift 
quasars, and could well be significant in the observed en- 
richment of the intergalactic medium of galaxy clusters. 
The fraction of metal to observed dust may thus be signif- 
icantly larger than the quoted value of 2.5. In conclusion, 
it is difficult to estimate a more precise value for Smd, 
which probably ranges between 1 and 10. 

If the radiative energy that is emitted along with the 
generation of heavy elements in a galaxy is released over 
a timescale rgp, the average total luminosity of the stars 
during that time would be 



6 X 10 



Em I TSF = 0.02 Smd Md c^ I rsF 

11 / Smd\ /I0^yr\ / 7.5 \ / 5250 \ ^ 



TSF 



Kd230 , 



) U mJy j 



(9) 



using Eqs. (6), (8) and (2). Would this energy be steadily 
emitted from the Big Bang to 2; « 4, i.e., over a time 
of ~ 1 Gyr, then the average luminosity would be simi- 
lar to that of local ULIRGs, which is bright, but below 
the detection limit of (sub)millimetre observations. It ap- 
pears unlikely that such a high rate of star formation has 
been maintained for a Hubble time (however, see Isaak et 
al. 2001). More likely, vigorous star formation started at 
a later time, in a series of short duration bursts, during 
which the mass of the heavy elements implied by the dust 
emission was produced. 

One-third of the PSS quasars studied in this paper 
have LpiR ~ 10^'^ Lq. If most of this luminosity comes 
from a starburst, Eq. (9) shows that its duration cannot 
much exceed r-^ 50 Myr. This is the typical duration of 
starbursts in local ULIRGs. The duration of the peak lu- 
minosity of the bright quasars is close to the Eddington- 
Salpeter time, ~ 50 Myr (e.g., Rees 1984). It is hkely that 
most of this peak-luminosity phase of the PSS quasars that 
we observed, took place in the studied redshift range 4-4.5 
(which spans a total time ^ 150 Myr), since few similar 
objects are known at z > 4.5. Because about one third 
of the PSS quasars were detected at 1.2 mm, the duration 
of their 10^^ L© starburst phase cannot be much shorter 
than 15 Myr. It could well exceed this value if similar star- 
bursts extended over the period when the quasar are less 
luminous, as possibly supported by the high detection rate 
of the low luminosity part of our sample (Fig. ||) . 

In summary, starbursts in high redshift quasars can 
consistently account for the far-IR luminosity, the ob- 
served dust content, and the fraction of quasars which 
are bright at far-infrared wavelengths. However, we can 
not exclude the possibility that the observed far-IR lu- 
minosity is due to AGN heating, and that the dust was 
generated in weaker (below our detection limit) and longer 
(to generate the required dust mass) past starburst events. 
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3.6. Star formation and black hole growth 

Massive black holes are now recognized as a normal, per- 
haps ubiquitous, component of elliptical galaxies and spi- 
ral galaxy bulges. There is increasing evidence that the 
growth of central supermassive black holes and the for- 
mation of the first stellar populations are closely linked 
(see, e.g., Combes 2001, and Haehnelt & Kauffmann 2000 
and references therein). Therefore the presence of strong 
starbursts in luminous high-redshift quasars would not be 
surprising. 

In nearby galaxies, there is good evidence for a corre- 
lation between the black hole mass, Mbh, and the bulge 
mass, Mbui, with 



Mbui/Mj 



BH 



160 - 500 



(10) 



(Magorrian ct al. 1998; Kormendy & Richstone 1995; 
Gebhardt et al. 2000 and references therein). Although 
the average value of this ratio is close to 500, it appears 
(see Laor 2001) that it is lower for the most massive local 
black holes. 

The z > 4 PSS quasars we observed are extremely 
luminous, with a typical blue magnitude Mb = —27.5 
(Table 1). Adopting the blue luminosity bolometric correc- 
tion factor of ~ 12 derived for the PG sample by Elvis et 
al. (1994), the typical bolometric luminosity of the quasars 
is Lboi ~ 10"^^ Lq. Only a fraction, e = O.leo.i, of the mass 
equivalent energy accreted onto the black hole is radiated, 
so that the rate of growth of a black hole relates to its lu- 
minosity 

70 Moyr-i Lboi 



Mbh = 



L 



bol 



e c^ 



(11) 



The ratio between the star formation rate and the black 
hole growth rate is thus 



Msf 
M, 



= 140 6mf <5sb eo.i (ipiR/iboi) 



(12) 



BH 



where we used Eqs. (4), (5) and (11). The FIR to bolomet- 
ric luminosity ratio of the millimetre- rfetecferf PSS quasars 
ranges between 0.04 and 0.2, with an average value of 0.1, 
whereas the average ratio over the whole sample is « 0.04. 
The rate at which stars form in the PSS quasars is there- 
fore more than a factor ten lower than expected if stars 
and black holes formed simultaneously with the same ratio 
of bulge to black hole mass, as observed in local galaxies. 
The fact that the far-IR luminosity in the sample does 
not appear to correlate with the blue luminosity (Fig. 2), 
indicates that the instantaneous rates of bulge star for- 
mation and of black hole growth are not directly related, 
and that the rate ratio of Eq. (12) may vary in time. 
However, there is a similar relation between the global 
masses Mbui/Afen- The mass of a black hole radiating 
near the Eddington limit is (e.g. Rees 1984) 



M 



BH 



2.6xlOHl0 (Lboi/lO^^Lo) 



(13) 



Since the bulge mass can be related to the mass of heavy 
elements released (Eq. 7), thus to the dust mass emitting 



the observed far-IR radiation (Eq. 8), which is propor- 
tional to the millimetre flux detected (Eq. 2) , we find that 



Mbui 
M, 



BH 



6(5mf 



'md 



10^ 
L 



bol 



7.5 



( 



5. 



250 



(14) 



^^230/ V™Jy 

which is independent of the time scale of star formation 
and of the source of dust heating. Despite the uncertainty 
on Smd and Kd230i this ratio is again an order of mag- 
nitude smaller on average than the local observed ratio 
Mbui/MsH ^ 160 (Eq. 10). 

The bulge to black hole mass ratio which is implied for 
the present sample of z > 4 PSS quasars thus supports the 
idea that supermassive black holes grow relatively faster 
at early epochs than their surrounding host galaxies. This 
has also been suggested by observations of host galaxies of 
Icnsed quasars at high z (Rix et al. 1999) and by modelling 
(Kauffmann & Haehnelt 2000). 

The relatively faster black-hole growth compared to 
the bulge growth at z > 4 could be in part due to the spe- 
cial conditions of the primordial gas infall into the central 
potential wells of very young galaxies and to the larger 
proportion of gas in mergings at z > 4 than at later epochs 
(e.g., Haehnelt & Rees 1993 and Rees 1998). 

Such effects could be enhanced in the PSS quasars 
which have bolometric luminosities much higher than 
those of the average quasar population. The high-redshift 
PSS quasars were selected by their optical brightness, so 
only the most luminous objects were discovered, which are 
not representative of the median z > 4 quasar population 
(Fan et al. 2001), nor, a fortiori, of the median z w — 2 
population (Boyle et al. 2000). The PSS quasars could 
thus be unusual objects which at later times have more 
massive black holes than typical quasars with compara- 
ble bulge mass, or they later accrete a larger stellar bulge 
mass without much growth of the black hole. 

3.7. Relation with other quasar properties 

There is as yet very little information published about the 
PSS quasars listed in Tables 1 and 2, in particular about 
their spectra, their environment and possible lensing. The 
information from the few published spectra (Kennefick et 
al. 1995a, b) is given in Tables 1 and 2. This information 
is too limited to allow us to confirm that broad absorp- 
tion lines (BAL) or/and broad and weak emission lines 
indicate strong millimetre emission (Omont et al. 1996a). 
Djorgovski (1998, 1999) have observed the surroundings of 
about twenty such quasars and found companion galaxies 
in almost every case. Specific information is only given for 
a few sources (as indicated in Tables 1 and 2). These find- 
ings give further support to the general arguments that 
quasars at 2 > 4 may form at the very highest peaks of 
the primordial density field, in the cores of future giant 
ellipticals (Djorgovski 1998). The presence of associated 
absorption systems points in the same direction, but we 
still lack conclusive information on this issue. Indications 
that the quasars detected at 1.2 mm have a rich environ- 
ment could be related to recent results on high-z radio 
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galaxies showing spatially distributed dust and CO emis- 
sion (Ivison et al. 2000b; Papadapoulos et al. 2000), to the 
discovery of clusters of Lya emission sources (Pentericci 
et al. 2000) and, if it is not due to lensing, to the dou- 
ble peak observed in the 1.2 mm dust continuum of BR 
B1202-0725 (Omont et al. 1996b). 

4. Conclusion 

The availability of ever more sensitive bolometer detec- 
tors and of vast optical surveys have lead to the discov- 
ery of thermal emission from tens of z > 4 radio-quiet 
quasars. (Sub)millimetre astronomy is thus entering an 
era of statistical analysis of star formation in high-z ultra- 
luminous quasars. The detection of warm dust with masses 
> 10^ M0 and far-infrared luminosities of ~ 10"'^^ Lq sug- 
gests substantial star formation activity that could be as- 
sociated with a large initial collapse or with major merg- 
ing. Millimetre and submillimetre observations thus pro- 
vide unique information with which to study the parallel 
growths of bulges and central black holes at very large 
redshift. The exceptional objects we discussed may trace 
the very highest peaks of the primordial density field, in 
the cores of proto-clusters and of future giant ellipticals. 
An important question still to be addressed is how the 
millimetre emission of bright quasars at z > 4 compares to 
that of quasars in the redshift range z ~ 1.5 — 3, at which 
the quasar activity peaks. In earlier studies (Omont et al. 
1996a; Guilloteau et al. 1999) only a few detections were 
achieved in this redshift range. 
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